clear all

%% report parameters that differ across models
% ordering: baseline, no slow recovery, crra, AR(1), time-varying s

model_file_names = {'run1a_baseline_sr.mat','run1b_baseline_no_sr.mat',...
    'run5a_crra_sr.mat','run7a_Nstate_sr_fit_v3.mat','run6a_time_varying_s_sr.mat'};

for i = 1:5
    load(model_file_names{i},'parms')
    model_parms{i} = parms;
end

%% check parameters that should be equal

% all models: iota

% all models except: 

% everything except AR(1) model: grid_z(1), StateTransitionProbs

%% report results

% productivity process: value in L regime
fprintf('\nproductivity process: value in L regime\n')
for i = 1:5
    if i~=4
        temp = model_parms{i}.grid_z(1);
    else
        temp = NaN;
    end
    fprintf('& %4.4f', temp)
end

% productivity process: sigma_z
fprintf('\nproductivity process: sigma_z\n')
for i = 1:5
    if i~=4
        temp = NaN;
    else
        temp = model_parms{i}.sigma_z;
    end
    fprintf('& %4.4f', temp)
end

% risk sharing: x_bar
fprintf('\nrisk sharing: x_bar\n')
for i = 1:5
    fprintf('& %4.3f', model_parms{i}.x_bar)
end

% risk sharing: sigma_x
fprintf('\nrisk sharing: sigma_x\n')
for i = 1:5
    fprintf('& %4.3f', model_parms{i}.sigma_x(1))
end

% time preference: beta
fprintf('\ntime preference: beta\n')
for i = 1:5
    fprintf('& %4.4f', model_parms{i}.beta)
end

% coeff of rra: gamma
fprintf('\ncoeff of rra: gamma\n')
for i = 1:5
    fprintf('& %4.2f', model_parms{i}.gamma)
end

% ES: consumption index
fprintf('\nES: consumption index\n')
for i = 1:5
    fprintf('& %4.4f', model_parms{i}.chi)
end

% unemployment benefits: b
fprintf('\nunemployment benefits: b\n')
for i = 1:5
    fprintf('& %4.4f', model_parms{i}.b)
end

% vacancy creation costs: kappa
fprintf('\nvacancy creation costs: kappa\n')
for i = 1:5
    fprintf('& %4.4f', model_parms{i}.kappa)
end

% worker's bargain power: eta
fprintf('\nworker''s bargain power: eta\n')
for i = 1:5
    fprintf('& %4.3f', model_parms{i}.eta)
end

fprintf('\n\n')